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Introduction 

Over the recent years, the studies of semileptonic charm 
decays have steadily been (re)gaining importance. Within 
the standard model (SM) the CKM matrix element mod- 
uli \Vcs\ and \ Vcd\ are currently best constrained indirectly, 
through CKM unitarity, and not by direct measurement. 
Recently however, deviations from unitarity predictions 
have been reported in leptonic decays of the Dg meson. 
Since these processes are helicity-supressed in the SM, the 
complementary semileptonic channels might offer the op- 
portunity to enlighten the issue. Like in the kaon and B- 
meson decays it is essential to reduce the theoretical error 
on the leptonic and semileptonic decays, in order to match 
the current and future experimental accuracy and therefore 
look for the potential discrepancies between the direct and 
indirect \ Vcq \ determinations. 

In addition, semileptonic charm decays offer a useful 
testing ground for theoretical tools. The bulk of the mea- 
sured experimental charm events is represented by non- 
leptonic Z?-decays, a consistent QCD-based description of 
which is still missing. Only a small fraction consists of 
leptonic and semileptonic decays which we know how to 
describe and -at least in principle- compute from the first 
principles of QCD. In this endeavor, numerical lattice QCD 
calculations need additional theoretical inputs to control 
systematic uncertainties and to be able to connect to the ex- 
perimental measurements. Conversely, calculations based 
on heavy quark expansion (HQE) and operator product ex- 
pansion (OPE) need to control power corrections, which 
are of utmost importance for example in the extraction of 
I Vub I from inclusive S-meson decays, and are expected to 
be far more pronounced in inclusive charm decays. 

Exclusive semileptonic decays 

Motivation 

Within the SM the I^s and Vcd CKM moduli can be de- 
termined from CKM unitarity yj 
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Figure 1 : The tree-level charged-current process contribut- 
ing to leptonic (left) and exclusive semileptonic (right) D 
decays 

lattice Clil one obtains 

\Vcd\L = |T4d|c/T(1.00±0.05), (2a) 
\Vcs\l = |ycs|c/T(1.08±0.03), (2b) 

in particular, there is a 2.3(7 tension between such deter- 
mination of I Vcs I and the CKM unitarity fit. This calls for 
cross-checks of lattice QCD calculations as well as experi- 
mental measurements. 

D Piu 

The differential decay width of the semileptonic decays 
of charmed (D, Dg) mesons to light pseudoscalar mesons 
(P — TT, K, r]^ ^) and light leptons (£ — fi, e) in the SM can 
be parametrized in terms of two kinematical variables 
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where — {p' — is the momentum exchanged 
squared, |p| — ^ X{mjj, q'^, mp)/2m£) is the absolute 
three-momentum of the final-state meson in the D rest- 



frame [A (a, b, c) = + 



— 2{ab + be + ac)] and 



9e is the angle between the directions of final state meson 
and lepton again in the D rest-frame. The non-perturbative 
QCD dynamics is encoded in the relevant quark current 
matrix element and can be parametrized in terms of the 
suitable form factors, e.g. 



\Vcd\uT = 0.22508 ±0.00082, (la) 
\Vcg\uT = 0.97347 ±0.00019. (lb) 

On the other hand, these quantities can also be extracted di- 
rectly from recent experimental measurements of leptonic 
decays ^ /^i^ L2J and Dg — ^ |3, 4, 5|. Using 

the averages of Cleo & Belle measurements taken fom |6| 
and the most precise decay constants calculations on the 
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{P{p)\qrc\D{p')) = Uiq' 



(4) 



where P = p' + p. In the SM the contribution of /o to the 
decay rate is (helicity) suppressed by the charged lepton 
mass and has thus been neglected in eq. Q. The relevant 
CKM matrix element | Vcq \ can be extracted from experi- 
mental measurement of the (partial) decay width provided 



(1) the form factor normalization is known at a single kine- 
matical point (customarily at = 0); (2) a (partial) phase 
space integral can be extracted and compared to experi- 
ment 191. The form factor normalization can be calculated 
e.g. using QCD (light-cone) sum rules ifTOl . the precision 
of which however is intrinsically limited. Therefore a more 
accurate determination can only possibly be expected from 
lattice QCD simulations. 

In the past, lattice calculations were only accurate 
in a kinematical region of diminishing phase-space near 
'irnax = ("^-D~"^p)'^- A subscqucut extrapolation between 
lattice and experiment was needed using suitable form fac- 
tor parameterizations. Recent lattice studies however can 
already cover the complete kinematical region in D de- 
cays ifTTl [T2I [T3I. Extrapolation procedures nevertheless 
remain important in exclusive determinations of | Vub \ from 
B ^ TT transitions with considerably larger phase space. It 
is illustrative to consider the present experimental precision 
of the measured B irtv partial branching fractions with 
a cut on close to where lattice QCD studies can provide 
the form factor normalization lfT4l 

BF{total) = 1.36 x 10-^(1 ± 0.05) , (5a) 
BF{q^ < 16GeV^) = 0.94 x 10"''(1 ± 0.07) , (5b) 
BF{q^ > 16GeV^) = 0.37 x 10^^(1 ± 0.10) . (5c) 

By only considering the high q^ region, the experimental 
precision is reduced by as much as 40%. In order to take 
advantage of the complete available experimental statistics, 
a careful control of theoretical errors, associated with form 
factor parameterizations is needed. These in term can be 
tested in the charm sector 

Exact (analytic) shape of the form factors is illusive, as 
it is non-perturbative in nature. Lattice QCD can be used 
to calculate the values of the form factors numerically at 
individual q^ points, leaving the question of how to best 
(extra-)interpolate between the various q^ regions open. 
Most form factor parameterizations employ an expansion 
in functions of q^ . 

• The simplest Taylor expansion around q^ = Q has 
proven useful in if ^ tt transitions, where the re- 
quired convergence radius is small g^a^. — {rriK — 
m^)^. It is however less suitable for semileptonic de- 
cays of heavy mesons, where the phase space is much 
larger. Perhaps even more important however is the 
presence of nearby resonance poles, close to the phys- 
ical region in D ^ P (and B ^ P) transitions. 

• An orthogonal approach is represented by the sum 
over t-channel resonance pole contributions of the 
form (g^ — mf)"^. It has been traditionally used both 
in theoretical form factor calculations as well as in the 
experimental fits. 

• Using a conformal mapping of the complex t plane 
and taking into account the position of the physi- 
cal cut at t+ = (mu + mp)"^, one can general- 
ize the expansion around any point in terms of 



— to). The required convergence radius is 
bounded by [0, |z„a.| < 1] (c.f. IB]). 

Some additional information is provided by the form factor 
dispersion relations. Namely, m D ^ K decays the first 
t-channel resonance (I?*) is below the physical cut at q^ = 
t+ (the same reasoning holds also for B ^ tt transitions 
with the B* resonance) 
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t — q'^ — ie' 
(6b) 



Close to the kinematical end-point at q^ax^ the /+ form 
factor is therefore expected to be completely dominated by 
the first resonance pole term (see also IIT6I ). 

Finally, in the limit nic 00, niq 0, the form factors 
obey the so-called heavy quark scaling relations. When the 
final state meson is soft in the D rest-frame (near q^ax)' 
heavy quark effective theory predicts IT?! 
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(7a) 
(7b) 



On the other end of the kinematical phase space (near q^ = 
0) similar scalling relations have been initially found by 
using QCD sum rules li 1 8i 
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(8) 



They were subsequently confirmed within the soft coUinear 
effective theory (SCET) |fT9ll20l which in this limit predicts 
also 
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(9) 



where Ep is the energy of the outgoing meson in the rest- 
frame of the initial D meson. These relations are broken by 
(potentially large) perturbative (mc), and power (A/mc) 
corrections. 

Phenomenologically, the single pole ansatz for f^{q'^) 
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(10) 



has been found inconsistent with measured spectra of D ^ 
■kIv and D Kiv already some time ago 11211 l22l . 

A simple three-parameter ansatz respecting leading- 
order HQ scaling relations was constructed by Becirevic 
& Kaidalov (BK) | 
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Considered as a truncated pole expansion, by construction 
it assumes 
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or equivalently a k, These assumptions were 

found not to hold for -> tt, L» ^ if decays II24II251 On 
the other hand, at present precision the data are still well 
represented by an unconstrained BK fit . 

Using dispersion relations and conformal t-plane map- 
ping, an alternative description of the form factor shape can 
be constructed ifTSl 
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where P{q^) subtracts away poles below the continuum 
threshold at t+ and (f>{q'^,to) is a normalization factor 
conventionally fixed from perturbative OPE unitarity (see 
also ||26|). Such parameterizations have been used exten- 
sively to describe semileptonic B decays. In D ^ P tran- 
sitions, existing data has already proven much more con- 
straining than unitarity bounds ||27J . Nonetheless, by com- 
paring to experiment while using the form factor nor- 
malization from the lattice lfT2ll one obtains for the relevant 
CKM moduli ffl 



|Vcd|sL 
\Vcs\sL 



|KrfkT(0.99±0.11), (14a) 
|K.|aT(1.05±0.13). (14b) 



At present precision (dominated by the theoretical errors on 
the lattice form factor normalization) the values obtained in 
this way are well consistent with CKM unitarity, although 
the central value of \Vcs \ is larger than one, as in the case of 
leptonic decays. Finally, once the theoretical lattice QCD 
error will approach the present experimental sensitivity, the 
consistency of the applied form factor parameterizations 
will also be put to the test. 

Beyond the SM, new operator contributions may con- 
tribute signicantly to the decay rate 
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(15a) 
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To describe possible tensor contributions, the knowledge of 
one additional form factor is needed (/t)- In both SCET & 
HQET limits [TLM however one finds 1.28. .29. ,30] 



(16) 



The relation can be traced back to heavy quark spin sym- 
metry relating the matrix elements of the a^,y and 7^ Dirac 
structures acting on a heavy quark field. However it also 
requires the HQ and SCET form factor scaling laws to 



hold. In addition, the fr form factor depends on the QCD 
anomalous dimension of the tensor operator (whereas the 
conserved vector current does not). Therefore one can 
only speculate that the relation should approximately hold 
throughout the kinematical region at some operator match- 
ing scale close to the charm mass. The proposal can of 
course be tested on the lattice. 

Phenomenologically, new operators may contribute dif- 
ferently to leptonic and semileptonic decays. In particular, 
while chiral current interactions produce similar enhance- 
ment in both leptonic and semileptonic modes (propor- 
tional to the /+ form factor in the semileptonic case), scalar 
or tensor interactions need to scale with the lepton mass in 
order to produce similar enhancement in r and n leptonic 
modes as indicated by experiments. This means that such 
interactions will contribute negligibly to the semileptonic 
decay rates (r channel is kinematically forbidden). Other 
observables, available in the semileptonic mode may help 
to discriminate among contributions. A particular asym- 
metry was constructed in lISTl 

T{Ei^ > 0) - T{Ei^ < 0) 



where 



r{Ee^_ > 0) + r{Ee^_ < 0) 
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(17) 



(18) 



and El is the energy of the outgoing charged lepton in the 
rest-frame of the initial D meson. This observable is di- 
rectly sensitive to the interference between SM and scalar 
or tensor interactions. A dedicated experimental analysis 
would possibly need to determine the feasibility of such a 
measurement. 

D Viu 

Traditionally the semileptonic decays of charmed 
mesons to light vector mesons (V = p, K* ,uj,(j)) were 
considered less interesting due to the larger number of the 
form factors needed to describe the decay rate, less theo- 
retical control over their evaluation, as well as more chal- 
lenging experimental analyses. On the other hand, the vari- 
ous polarizations of the vector meson in the final state give 
access to more observables suitable for distinguishing be- 
tween SM and possible new physics contributions. 

In the SM (as well as in presence of new scalar contri- 
butions) the D ^ V transitions can be fully described in 
terms of four form factors V, ^0,1.2 defined through 
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If additional tensor contributions are present, these require 
the knowledge of additional three form factors Ti_2,3 
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While lattice QCD computations can provide normaliza- 
tion of the form factors at various kinematical points, HQ 
scaling laws and relations can again be used together with 
dispersion relations to construct useful parameterizations. 
In the soft V limit lEl 
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(21a) 
(21b) 
(21c) 
(21d) 
(21e) 



while near q 
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all form factors should scale as 
rrijj ' lfT9l l20l . There are also additional form factor rela- 
tions, broken by q;s(toc), A/toc corrections, which allow 
analogues of the BK parameterization for H ^ V decays 
to be constructed I 
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(22b) 
(22c) 
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with, y(0)/Ai(0)/r(0) w l,aW « - 0(1). The 
different nearest resonance poles in V, Ai and Ti are re- 
flected in the different normalization of the variables x,y — 
Q^/'i^poie- Note that eqs. ( I22ab and ( I22el ) already tacitly 
hint at the relation 

Ti^V, (23) 



throughout the allowed kinematic region, which can be 
traced back to the same origins as eq. ( fTSI ) but is expected 
to receive additional large corrections due to non-negligible 
light vector masses. The differential decay width is usually 
written in terms of Helicity amplitudes o 
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which can in term be related to certain combinations of the 
various D ^ V form factors. Experimental information 

on the shapes of o is already available |f34l [35] [36l 

and can be used to test form factor parameterizations like 
the z-expansion |24| or the modified (BK) pole ansatz fJTll . 
In addition to the (partially) integrated decay rate, suitably 
constructed asymmetries, analogous to (Til might provide 
more direct access to anomalous scalar, tensor contribu- 
tions. 

Inclusive Semileptonic Decays 

Introduction 

Over the last years, there has been a tremendous progress 
in the determination of the | Vub \ CKM matrix element from 
the measurement of the inclusive decay rate B XuiP- 
Using OPE and HQE the theoretical predictions ||40l |4T] 
|42l|43l|44l have reached a precision below 10% (Wl. Es- 
timation of relevant power-supressed non-perturbative op- 
erator matrix element values from inclusive analyses of 
B Xctv has been instrumental in this effort. At the third 
order in the heavy quark mass expansion {l/m\) however, 
effects due to dimension-6 four quark operators, the so- 
called weak annihilation (WA) contributions appear, which 
cannot be extracted from the inclusive 6 — > c analysis. In 
addition, these are IGtt^ phase space enhanced compared 
to LO & NLO contributions (such enhancement does does 
not seem to appear at dimension-7 [45 1). 

Recently, the inclusive semileptonic decay rates of 
charmed mesons {D XlD ) have been determined ex- 
perimentally, yielding for the branching fractions ISS] |39l 

B{D+ ^ Xeu) = (16.13 ±0.20 ±0.33)%, (25a) 
B{D°^XeP) = (6.46 ±0.17 ±0.13)%, (25b) 

with similar results for muons. On the theory side, one 
may attempt to accommodate these numbers by treating the 
charm quark mass as heavy, and perform an OPE calcula- 
tion by expanding in as{mc), A/rric- In the process, one 
needs to estimate contributions of local operator matrix ele- 
ments, which can possibly be related to the ones appearing 
in inclusive B decay analyses or estimated on the lattice. 

Heavy quark and operator product expansion 

The total inclusive decay rate can be related to the ab- 
sorptive part of the forward scattering amplitude by apply- 
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Figure 2: Examples of contributions to the transition oper- 
ator T (left) and to the corresponding local operator (right). 
The open circles represent the insertions of the weak effec- 
tive Hamiltonian. The full circles represent the insertion of 
a local AQ = operator. 



ing the optical theorem 
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where 



T = Im z / d^xT[Heff{x)Heff{0)] 



(26) 



(27) 



and T is the time-ordering operator while Hef / is the effec- 
tive weak Hamiltonian mediating charged current semilep- 
tonic processes. For a decaying hadron containing a sin- 
gle heavy quark one can, assuming quark-hadron duality, 
evaluate the quantity by expanding it in powers of the in- 
verse heavy quark mass, systematically including pertur- 
bative QCD corrections. Applying the strategy to semilep- 
tonic D meson decays one obtains Il46ll47l 

r(D-«i;) = 2M|„„|»g,,) 
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where r — rnl/m^, gg is the QCD coupling constant, 
G^iy — [DjT)], where T> is the QCD covariant derivative, 
g{r) is the leading kinematical factor, while /o,i (''') contain 
relative kinematical as well as perturbative as{mc) cor- 
rections. The ellipses denote additional 1/m^ and higher 
power corrections which are not phase-space enhanced. 

The non-perturbative operator matrix elements appear- 
ing in eq. ( l28T l can be parameterized using the heavy quark 
equations of motion 
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where 
eters 



V ~ v{v ■ V), and related to the HQE param- 
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Similarly, the operator matrix elements appearing at 1/ m^ 
can be parameterized in terms of two additional parameters 
(Pd.sl) ll48l l49l . The /i^r and parameters depend on 
the employed heavy quark mass scheme and the associated 
matching scale |48|. Finally, also at 1/m^ the contributions 
involving light flavors 



Q7M(i-75)g'gV(i-75)Q, 
Q(i-75)9Y(i-75)Q. 



(32) 



are conventionally parametrized in terms of deviations 
from the complete factorization or vacuum saturation ap- 
proximation by introducing the suitable bag parameters 
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In this way one defines 2mp)By^A — {D\ Ov-a — 
Os-p \D) = f^rrijj{Bi — B2). These quantities depend 
on the scale (/^wa) at which the matrix elements are eval- 
uated and can mix with phase space non-enhanced contri- 
butions (pn) ll50l . The early estimates of B12 were done 
in the framework of QCD sum rules fsT, '52^ ^ 53l . They 
have been also computed on the lattice |54 5 51 l56l l57ll . 
However, since possible contractions of the Ov-a.s-p 
operators involve closed fermionic line topologies (eye- 
contractions) ||58]| . such lattice estimates are necessarily in- 
complete I 
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Figure 3: Graph topologies of WA operators contributing 
to inclusive decays of heavy hadrons. The full circles rep- 
resent the insertion of a local AQ = operator Only the 
left topology is estimated on the lattice. 

Putting all the pieces together one obtains for the decay 
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where the leading /iwA scale dependence (which should 
cancel between the two terms in the fourth line) is 
given by d(r,^wA) [-lOr^ + 32r^ - 2Ar'^ - 32r + 
241og (/^wa/'^c) + 34]/3. The expressions for g{r) and 
Pi,2(?') can be read from the Xi functions in ref. ||59l as 
g(r) =Xo,pi,2(r) -4X1^2/3X0. 

Phenomenological analysis 

The most sensitive parameter entering the inclusive de- 
cay rate analysis is the charm quark mass. Without impos- 
ing any kinematical cuts it enters the total decay rate with 
the fifth power. From lattice II6OI , Charmonium SR 11671 
l62i . b ^ c spectral fits ll63l one obtains a fairly accurate 
and consistent {MS) value of mc(mc) = 1.27(2)GeV. 
Values of and pd,sl can also be extracted experi- 
mentally. Namely /i^ = (3/4) [m|,. - m|,] = 0.41GeV^, 
while w 0.4GeV^ pi, « 0.2GeV^ and p^^ « 
-0.2GeV^ are obtained at the 10% - 20% precision from 
a fit to the b ^ c spectrum |64| in the so-called kinetic 
heavy mass scheme [651 at the scale of 1 GeV. Adding 
up all these known contributions in eq. ( |34] |. the experi- 
mental values are saturated to sa 70%. The calculation ex- 
hibits a very slow perturbative & power convergence, the 
details will be presented elsewhere |66|. In light of sizable 
residual scale and scheme dependencies (the best perturba- 
tive convergence is achieved for a very low kinetic charm 
mass scale of 0.5 GeV) the result should thus be taken as 
tentative at best. Assuming the unknown enhanced 1 / 
terms saturate the rate, one can in any case obtain a guesti- 
mate |67| of (using fo = 208(4)MeV Q) 



{Bi - B2)(2GeV) w 0.05, 



(35) 



which appears to be a reasonably small number One 
should note however, that a fit to the inclusive _D+ " de- 
cays can only probe spectator see-quark operator contri- 
butions (where the flavour q' appearing in eq. (l32T i does 
not match the (light) flavour of the initial charmed hadron), 
while the valence quark contribution is only relevant in 
Ds decays. Once determined, this could be related to 
the B+ XJ+v and B" XJ+v width difference 
via SU{3) symmetry f58 |. For more accuracy and cross- 
checks, improved lattice QCD estimates of Bi^2 are of 
course called for. 



In the meantime the phenomenological analysis could 
certainly also be improved by considering higher spectral 
moments in addition to the total decay rates (as done in 
b ^ c case) in order to leverage more control over the 
OPE convergence. Cleo already published the lepton mo- 
mentum spectra in inclusive semileptonic D decays with 
a lower cut [38| which however requires a more involved 
treatment. Utilizing the moments, a more direct access 
to power corrections and even possible duality violations 
(c.f. iMlES) could be obtained, since e.g. WA contribu- 
tions are expected to dominate near the spectrum end-point, 
c.f. iCOllTll. 



Conclusions 

The exclusive D PiD decays offer the opportunity 
to confront the recent puzzles concerning the determina- 
tion of \Vcq\ CKM matrix elements from leptonic charm 
decays. Here the most critical input from theory is to pro- 
vide an accurate normalization of the relevant form fac- 
tors. While new results from lattice QCD calculations 
are eagerly awaited, presently the extraction of the afore 
mentioned CKM elements is not yet competitive with lep- 
tonic decays. Furthermore, various form factor parame- 
terizations used to model semileptonic D decay spectra 
are also employed in B decays. By comparing the mea- 
sured D decay spectra with lattice results and fitting both 
to these parameterizations one can test extrapolation pro- 
cedures needed in exclusive determination of \Vub\- Fi- 
nally, by using even more of the present experimental in- 
formation available, anomalous contributions due to scalar, 
right-handed or tensor currents could be probed. Impor- 
tantly, such an analysis introduces almost no new hadronic 
uncertainties. In this respect the D — > V£D decays offer 
potentially even more useful observables, sensitive to NP 
contributions, with the caveat of heaving to estimate more 
hadronic form factors. 

Inclusive D Xlv decays offer the possibility to test 
and extract power suppressed spectator contributions to in- 
clusive semileptonic decay widths of heavy hadrons from 
experiment. These are relevant for a reliable extraction of 
\Vub\ from inclusive B decay measurements. Existing es- 
timates would need to be confronted with improved lattice 
QCD calculations, while the convergence and validity of 
the OPE in the charm sector needs to be examined care- 
fully. More experimental observables are already available 
than the total rate and dedicated experimental and theoreti- 
cal analyses of these are called for. 
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